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Abstract 

We propose an alternative approach that avoids the nonlinear 



' equations for the Fourier coefficients that appear in the method of 

o" 

harmonic balance. We apply it to two simple illustrative examples. 



1 Introduction 

Since one-dimensional nonlinear oscillators are treated as simple pedagog- 
ical examples in many standard textbooks [1-5] one is tempted to believe 
that they are trivial problems. However, from time to time there appears to 
be a renewed interest in them as shown by a recent unbelievable article [6]. 
The equations of motion for the simplest nonlinear oscillators are suitable for 



illustrating the application of several approximate methods that include vari- 
ous forms of perturbation theory, variational approaches, Fourier expansions, 
etc. [1-5]. 

The method of harmonic balance [2], which is based on the Fourier ex- 
pansion of the solution of a periodic nonlinear oscillator, yields accurate 
results but requires the solution to a rather complicated system of nonlinear 
equations for the expansion coefficients. It is not that difficult to select the 
physical solution of a nonlinear equation in one variable, but the problem 
becomes increasingly complicated as the number of equations and variables 
increases (specially if there are many unwanted solutions besides the physical 
one). The purpose of this communication is to propose a simpler approach 
that avoids the nonlinear system of equations for the Fourier coefficients. It 
is based on a systematic generalization of the approach proposed by Ren and 
He [6] in a paper plagued by errors and misconceptions [7]. In Sec. [2] we 
develop the method and in Sec. [3] we apply it to two simple but nontriv- 
ial illustrative examples. Finally, in Sec. H] we discuss the advantages and 
limitations of the approach. 

2 The method 

Consider the dimensionless equation of motion for a nonlinear oscillator 

u"{t) + f[u(t)] = 0, u(0) = A, «'(0) = (1) 

where u(t) is the displacement from equilibrium and the prime indicates 
differentiation with respect to the independent variable t. For simplicity we 
assume that the force —f(u) is an odd function of the displacement: f(—u) = 

2 



—f(u) and that the motion is bounded with a period T = 2tt/uj, where 
u is the angular frequency. The potential-energy function V(u) (f(u) = 
dV(u)/du) is, consequently, an even function of the displacement (V(—u) = 
V(u)). Such a symmetry and the initial conditions determine that w(T/4) = 
and therefore we can restrict our calculations to the interval < t < T/4 
where A > u(t) > 0. Obviously, if we know u(t) in this interval then we 
know it for all values of t. 

If we try to solve the equation of motion ([T]) by means of a truncated 
Fourier-series expansion 

TV 

u W(t) = Y^Aj cos[(2j - l)ut] (2) 
j=i 

and the method of harmonic balance [2] then we are led to a system of 
nonlinear equations for the expansion coefficients Aj. In what follows we 
propose a simple strategy that avoids such a problem. Notice that the ansatz 
(ED satisfies uM(T/4) = 0. 

If we define u^ 2n \i) = d 2n u{t) / dt 2n , u^\t) = u(t) then we can write 
the set of equations -u^ 2n )(0) = G n (A), n = 0,1,..., where G (A) = A, 
G X {A) = -f(A), G 2 (A) = f(A)f'(A), etc. Substitution of Eq. © into the 
left-hand side of these equations for n — 0, 1, . . . , N — 1 leads to the linear 
system of equations 

JV 

J2(2j - l) 2n u 2n A j = G n (A), n = 0, 1, . . . , JV - 1 (3) 

i=i 

that gives us the coefficients Aj in terms of the amplitude A. The matrix of 
this system of equations is the transpose of a Vandermonde one and therefore 
its determinant is nonzero. 
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If we now integrate the equation of motion ([T]) twice and take into account 
the initial conditions we have 

u(t)=A- f [ f{u{t"))dt" dt' (4) 
Jo Jo 

that we can apply iteratively. Here we carry out only one iteration 

v,W (t)=A- f f f(u [N] (*")) dt" dt' (5) 
Jo Jo 

and obtain the period from the resulting approximate solution u^(t) and 
the condition derived above: 

u [N \T/A) = tt (6) 

For N = 1 we simply have itW(t) = Ai cos(ut) and A\ = A from the only 
equation in ([3]). The system of equations with N = 2 yields 

9Auj 2 - f(A) 



A x 



lu 2 



A - f(A) ~ (7) 

These two approaches of first (N = 1) and second (N = 2) order are suffi- 
ciently accurate as suggested by the examples below. 

If the force is not an odd function of the displacement then we choose the 
ansatz u^(t) = A + A\ cos(u;£) + A 2 cos(2cut) + . . . + A N cos(Nujt), which 
should satisfy u^(0) = A and v\ N \T /2) = B, and resort to the independent 
condition v\ N \T/2) = B, where B is the other turning point. 
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3 Examples 



Our first example is one of the most widely studied nonlinear oscillators: the 
Duffing equation given by the function [1-5] 

f(u) = u + eu 3 (8) 

It is not difficult to prove that p = eA 2 is the relevant parameter of the 
model. That is to say, u(t)/A and other properties of the system, like, for 
example, the period, depend on A and e only through p and not separately 
or in other combinations. 

In the first approximation we derive the period [6] 

^M-T&H (») 

In order to verify its accuracy for large values of p we recall that the exact 
period satisfies 

lim y/pT(p) =T QO = 7.416298709 (10) 

p— >oo 

The approximate expression (jHJ) yields the reasonable estimate 

lim v ^rW(p)=TW« 7.12 (11) 

p^oo 

For small values of p we have 

T(p) = 2n- 3 ^ + 5 ^ + ... (12) 
\H) 4 128 K ' 



and 



T W = 2n - + . . . (13) 



We appreciate that (p) does only give us the leading term of the small-p 
series. 



For N = 2 we easily obtain 

125T 8 p (p + l) 3 - 7656vr 2 T 6 p (p + l) 2 + 120vr 4 T 4 (p + 1) (1607p + 735) 
-64vr 6 T 2 (48851p + 55125) + 12700800tt 8 = (14) 

The approximate period is given by the only real root T^(p) that satisfies 
T^(0) = 2n. It is clear that we have not avoided the nonlinear nature of the 
problem completely but we have reduced it to a an equation with just one 
variable which simplifies the task of choosing the physical solution. 

Without difficulty we derive the following expression for the approxima- 
tion of second order to 

125T^ - 76567r 2 T^ + 192840tt 4 T ( 4 2 - 3126464tt 6 T ( 2 + 12700800tt 8 = (15) 

from which we obtain J 1 ] 2 , 7.44 that is in fact more accurate than T^}. 
The Taylor expansion of Tt 2 l(p) about p = 

r 9l . , 37T/0 14317T , . 

T [2] p =2tt - + + ... 16 

KHJ 4 3200 V ; 

exhibits the first two terms of the exact small-p series. 

Fig. [1] shows the exact period and the two approximations just derived 
for small and moderate values of p. We appreciate that T^(p) is more 
accurate than (p) for all values of p as expected from the analytic results 
derived above. It is encouraging that the accuracy of the method increases 
notoriously with the order and that just the lowest two approximations are 
in good agreement with the exact result. 

Our second example is given by 

f{u) = u + e\u\u (17) 



and we realize that the relevant parameter is p = eA. This problem is as 
simple as the preceding one and we easily obtain = 6.869261369. 

A straightforward calculation shows that the approximate method out- 
lined in Sec. [2] yields 



T[1(P) v/p(vr 2 + 4) + 16 (18) 

and T£ ] fts 6.75 is a reasonable estimate of the exact large-p limit. These 
results were obtained by Ren and He [6]. 
For the small-p series we have 

T(p)=2n-^- + ... (19) 



and 



/ \ tt(tt 2 + 4)p , , 

16 



For the second approximation we have 

T 6 p(p + 1) 2 (9tt 2 - 16) - 8tt 2 T 4 (p + 1)[p(45tt 2 - 16) + 256] 

+ 167r 4 T 2 [p(3697r 2 + 1136) + 5120] - 294912tt 6 = (21) 

and the expression for the large-p limit 

r^(97T 2 -16)+87r 2 T ( 4 3 (16-457r 2 )+167r 4 r 2 (3697r 2 +1136)-2949127r 6 = (22) 

The estimate t£ ~ 6.866 is more accurate than T^j. 
The small-p expansion is also more accurate 

m rai/ x 7r(97r 2 + 20)p , , 

T (P) =2n ~ 128 + • • • (23) 

but neither T' 1 ' nor T' 2 ' gives more than the leading term exactly. 



Fig. [2] shows the exact period and the approximations of first and second 
order for small and moderate values of p. Once again we appreciate that the 
accuracy of the method increases with the order. 

4 Conclusions 

We have developed a straightforward systematic approach for the calculation 
of the period of simple nonlinear oscillators. It is particularly useful if we 
can integrate the equation ([5]) analytically which we easily do when f(u) is 
a polynomial or when a few terms of the Taylor series about u = give a 
reasonable approximation. 

Although present approach becomes increasingly complicated with the 
order, we think that its equations are always considerably simpler than those 
coming from the method of harmonic balance [2] . 
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Figure 1: Period for the Duffing oscillator. Circles, dashed and solid lines 
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Figure 2: Period for the discontinuous oscillator. Circles, dashed and solid 
lines stand for exact, first and second approximations, respectively 
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